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SUMMARY 


One  of  the  most  attractive  modifications  for  fire  tube  boilers  is 
the  introduction  of  "turbulators"  into  the  boiler  tubes.  These  inserts 
increase  gas-side  heat  transfer  coefficients,  thereby  increasing  the 
heat  removed  from  the  flue  gases.  Assuming  that  the  equipment  is 
properly  operated  and  adjusted,  the  boiler  efficiency  will  be  higher. 
Although  several  types  of  twisted  metal  strips  are  quite  widely  used  as 
original  equipment  or  as  retrofit  hardware,  heat  transfer  coefficients 
and  friction  factors  are  not  available  in  the  open  literature.  This 
information  was  obtained  in  the  present  study  for  three  popular  turbulator 
inserts . 

An  electrically  heated  flow  facility  was  developed  to  deliver  hot  air 
to  a  water-cooled  steel  tube  of  2.675  in.  i.d.  and  71.75  in.  length.  Tube 
wall  temperatures,  fluid  bulk  temperatures,  and  flow  rates  were  measured. 
The  cooling  coil  was  segmented  so  that  sectional  average  heat  transfer 
coefficients  could  then  be  derived  for  four  regions  of  the  tube. 

Reference  data  for  the  empty  tube  are  in  excellent  agreement  with 
the  usual  correlations.  Two  commercial  turbulators,  consisting  of 
narrow,  thin  metal  strips  bent  and  twisted  in  zig-zag  fashion  to  allow 
periodic  contact  with  the  tube  wall,  displayed  125  and  157%  increases 
in  heat  transfer  coefficients  at  a  Reynolds  number  of  10,700.  Insert 
orientation  had  at  most  a  10%  effect  on  the  heat  transfer  coefficient. 

A  third  commercial  turbulator,  consisting  of  a  twisted  strip  with  width 
slightly  less  than  tube  diameter,  provided  a  60%  increase  in  heat  trans¬ 
fer  coefficient.  The  friction  factor  increases  accompanying  these  heat 


transfer  coefficient  increases  were  1110%,  1000%,  and  120%,  respectively. 
These  data  should  be  helpful  in  assessing  the  overall  performance  gains 
to  be  expected  when  the  turbulators  are  used  in  actual  boilers. 
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INTRODUCTION 


Background 

The  U.S.  Army  is  activ  iy  pursuing  hardware  and  strategies  which 
will  lead  to  the  reduction  of  tael  oil  used  for  space  heating.  Even¬ 
tually,  many  of  the  oil-fired  plants  will  be  replaced  with  systems 
fired  by  coal.  In  the  interim,  however,  it  is  important  to  minimize  oil 
consumption . 

One  of  the  most  attractive  boiler  modifications  is  the  introduction 
of  "turbulators"  into  the  boiler  tubes.  One  popular  configuration  is 
a  narrow,  thin  metal  strip  bent  and  twisted  in  zig-zag  fashion  to  allow 
periodic  contact  with  the  tube  wall.  Twisted  ribbons  or  tapes,  of  tube 
diameter  width,  have  been  used  for  a  century.  Another  simple  insert  is 
a  coiled  wire  similar  to  a  spring.  These  inserts  increase  the  gas-side 
heat  transfer  coefficients,  thereby  increasing  the  heat  recovery  from 
the  flue  gases.  Assuming  that  the  equipment  is  properly  operated  and 
adjusted,  there  will  be  a  higher  boiler  efficiency.  Of  course,  boiler 
performance  will  reflect  the  increased  flow  friction  of  the  turbulators. 

Turbulator  inserts  are  advertised  in  rather  glowing  terms,  with 
typical  claims  of  up  to  15%  savings  in  fuel  based  on  field  operational 
data.  Unfortunately,  the  accuracy  of  such  field  tests  is  open  to 
question.  Furthermore,  it  does  not  appear  that  heat  transfer  and  fluid 
flow  characteristics  of  the  major  turbulator  inserts  have  been  reported 
in  the  open  literature.  For  these  reasons  and  the  fact  that  boilers 
for  military  installations  operate  at  low  heat  fluxes,  the  Army, 
through  the  U.S.  Army  Facilities  Engineering  Support  Agency  (FESA)  has 
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embarked  on  an  evaluation  program  to  determine  the  actual  energy  savings 
to  be  expected  from  the  installation  of  turbulators. 

It  is  the  intent  of  this  evaluation  program  to  estimate  what  the 
fuel  reduction  "could  be"  for  any  boiler  on  a  seasonal  basis,  rather 
than  what  the  reduction  is  during  a  short  period  of  steady-state  oper¬ 
ation.  Rather  than  test  an  actual  boiler,  the  Army  program  will  use 
National  Bureau  of  Standards  computer  models  to  simulate  boiler  per¬ 
formance.  As  preliminary  input,  typical  published  values  of  the  gas- 
side  heat  transfer  coefficient  and  friction  factor  for  twisted  tapes 
and  wire  coils  were  utilized.  Actual  values  of  the  thermal-hydraulic 
performance  of  commercial  turbulators  were  then  measured  in  the  Iowa 
State  University  Heat  Transfer  Laboratory.  The  measurements  were  car¬ 
ried  out  for  a  smooth  tube  and  for  the  same  tube  with  various  inserts 
under  conditions  of  cooling  of  a  hot  gas  at  flow  rates  representative 
of  actual  boilers.  This  report  concerns  the  findings  of  the  Iowa  State 
Heat  Transfer  Laboratory  experiments  for  three  commercially  available 
inserts . 


Turbulators  to  be  Tested 

Three  turbulator  inserts  were  provided  to  the  Heat  Transfer  Lab¬ 
oratory  by  FESA.  These  inserts  are  offered  commercially;  however, 
detailed  thermal-hydraulic  performance  data  do  not  appear  to  be  avail¬ 
able  in  the  open  literature.  The  manufacturers  or  distributors  of  the 
inserts  and  their  addresses  are  given  in  Appendix  A. 
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The  somewhat  complex  geometric  characteristics  of  the  inserts  are 
given  in  Figs.  1  through  3,  where  the  nomenclature  used  is  noted  on 
each  figure  and  a  letter  symbol  is  used  to  identify  each  insert. 
Numerical  values  of  the  various  geometric  dimensions  are  given  in 
Table  I.  Insert  pitch  is  defined  in  various  ways  in  the  references 
cited  in  this  report.  In  order  to  have  a  uniform  basis  for  comparison, 
this  report  uses  diametral  pitch:  the  number  of  inside  tube  diameters 
along  the  tube  axis  in  which  one  full  twist  (360°)  of  a  twisted  tape 
or  one  full  cycle  of  a  bent  strip  turbulator  occurs. 


Table  I.  Insert  characteristics. 


Property  _ Insert 


Type  A 

Type  B 

Type  C 

Fig.  1 

Fig.  2 

Fig.  3 

Length 

in . 

72.5 

73.0 

74.5 

Tube  span 

in . 

2.62 

2.56 

2.60 

Pitch 

(for  2.675  in. 
i.d.  tube) 

4.02 

2.80 

10.47 

Width 

in . 

0.880 

0.755 

2.60 

Thickness 

in. 

0.075 

0.085 

0.055 

Longi tudina 1 

Angle 

deg 

122° 

100° 

- 

Primary  Twist 

Angle 

deg 

±25° 

30° 

Secondary  Twist 
Angle 


deg 


-30° 


LITERATURE  REVIEW 


A  literature  search  on  the  effects  of  tube  inserts  of  various 
types  suitable  for  augmentation  of  gas  heat  transfer  coefficients  in 
fire  tube  boilers  was  conducted  using  the  current  Computerized  Bibliog¬ 
raphy  on  Augmentation  of  Convective  Heat  and  Mass  Transfer  [1].  Inserts 
most  applicable  to  the  present  investigation  of  fire  tube  boiler  turbu- 
lators  are  forms  of  twisted  tapes  and  wire  coils.  Some  early  studies 
report  complete  boiler  tests  with  and  without  twisted-tape  turbulators. 
The  incentive  was  to  confirm  on  a  stationary  boiler  the  advantages  of 
the  marine  practice  of  using  "retarders." 

Whitham  [2]  in  1896  inserted  "loosely  fitting"  and  gently  twisted 
(1  turn  in  10  ft)  tapes  in  the  44-4  in.  horizontal  tubes  of  a  100  hp 
boiler.  The  saving  of  fuel  with  turbulators  was  as  much  as  8%,  with 
the  improvements  being  observed  only  when  "the  boiler  plant  is  pushed 
and  the  draught  is  strong."  These  general  conclusions  were  reinforced 
over  50  years  later  by  Kirov  [3]  who  fitted  an  Economic  Boiler  having 
76-2.74  in.  i.d.  tubes  with  2.5  in.  width  twisted  tapes  of  1  turn  in 

1.17  ft  or  1  turn  in  2.33  ft.  With  full-length  tapes,  the  increase 

1  2 

in  boiler  efficiency  was  5  to  7  percentage  points  ,  which  was  equi¬ 
valent  to  a  fuel  saving  of  7  to  10%.  The  increase  in  fan  horsepower 
was  negligible  for  the  forced  draft  system  under  test;  however,  use  of 
the  retarders  in  a  natural  draft  system  was  not  advised. 

^Boiler  efficiency  is  defined  as  the  heat  transferred  to  the  water 
divided  by  the  heat  available  from  the  fuel  consumed. 

2 

Based  on  higher  heat  of  combustion  of  the  coal,  as  fired. 


In  a  recent  study  conducted  for  DOE  at  Brookhaven  National  Lab¬ 
oratory,  it  was  found  that  the  addition  of  turbulators  to  residential 
oil-fired  boilers  resulted  in  fuel  savings  of  2-8%  [4].  Details  are 
not  given  in  the  report. 

Actual  heat  transfer  and  flow  friction  data  found  for  twisted 
tapes  and  wire  coils  in  air  flows  are  described  below;  these  are  the 
only  cases  found  with  air  as  the  working  fluid  and  with  cooling.  Other 
data  using  water  as  the  working  fluid  have  been  obtained,  but  they  are 
not  considered  useful  to  this  investigation. 

Wire  Coils 

Early  work  on  this  type  of  enhancement  was  performed  by  Colburn 
and  King  (5|  on  conically  shaped  copper  wire  coils.  Wire  was  wound  in 
a  three-dimensional  conical  helix,  beginning  with  a  coil  of  essentially 
the  same  circumference  as  the  inside  of  the  2.625  in.  i.d.test  section 
and  progressing  to  the  apex  of  the  cone  in  ever-decreasing  radii. 

The  altitude  of  the  cone  thus  formed  was  about  4  in.  Nine  of  these 
coils  were  placed  in  the  3  ft  long  test  section  in  an  apex-to-base 
sequence.  No  information  is  given  on  wire  diameter,  nor  on  the  direc¬ 
tion  of  flow  relative  to  the  cone  sequence  in  the  test  section. 

The  test  fluid  was  air,  heated  by  electrical  resistance  heaters. 
The  test  section  was  water-cooled  by  a  helically-wound  copper  tube 
soldered  to  the  outside,  with  end  losses  and  axial  conduction  losses 
reduced  by  using  guard  coolers  on  sections  of  tube  just  upstream  and 
downstream  of  the  test  section.  Tube  wall  temperatures  were  measured 
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on  the  outside  wall  of  the  tube  at  five  equally-spaced  axial  locations 
by  the  thermocouples  embedded  in  the  solder  between  copper  cooling 
coils.  Entering  air  temperature  was  measured  by  obtaining  an  average 
of  temperatures  read  from  a  thermocouple  traversed  along  a  diameter  of 
the  tube.  The  exit  gas  temperature  was  taken  from  thermocouple  read¬ 
ings  in  an  insulated  mixing  chamber  at  the  test  section  outlet.  Inlet 
air  temperatures  ranged  from  107  to  390°C,  and  inlet  wall  temperatures 
varied  from  9.2  to  11.5°C. 

Data  are  reported  as  curves  of  h^/c^  and  pressure  drop  corrected 

for  kinetic  energy  losses  as  functions  of  mass  velocity  through  the 

test  section.  The  data  are  also  given  in  tabular  form.  Recast  in  the 

0  4 

more  familiar  form  of  NUjj/(Pr)  ‘  and  f  as  functions  of  Re^  based  on 
average  fluid  properties,  the  data  are  shown  in  Figs.  4  and  5. 

Twisted  Tapes 

An  early,  but  comprehensive  study  of  twisted-tape  inserts  was 
reported  by  Royds  (6j.  Electrically  heated  air  was  allowed  to  flow¬ 
through  a  tube  of  2.31  in.  i.d.  and  6.58  ft  length.  The  last  section 
was  cooled  by  water  flowing  in  a  shell  around  the  tube.  The  inlet  air 
temperature  was  measured  by  a  calibrated  platinum  resistance  thermometer. 
Other  temperatures  were  obtained  by  calibrated  mercury-in-glass  thermom¬ 
eters.  A  large  number  of  tests  were  made  with  inlet  air  temperatures 
varying  from  430°F  to  1000°F. 

Nine  different  sooted  twisted  tapes  of  1.94-in.  width  were  tested 
with  pitches  varying  from  4.72  to  38.5.  The  results  obtained  show 


1  WIRE  COIL  (COLBURN 

2  TWISTED  TAPE  PITCH 

3  TWISTED  TAPE  PITCH 

4  TWISTED  TAPE  PITCH 

5  TWISTED  TAPE  PITCH 

6  TWISTED  TAPE  PITCH 

7  TWISTED  TAPE  PITCH 

8  SMOOTH  TUBE 

_ I _ I _ l_ 

3  4  5 


AND  KING)  [6] 

*  2.67  (COLBURN  AND  KING)  [6] 
-2.33  (EVANS  AND  SARJANT)  [7] 
-4.67  (EVANS  AND  SARJANT)  [7] 
-4.72  (ROYDS)  [5] 

-  38.5  (ROYDS)  [5] 

-3.15  (THORSEN  AND  LANDIS)  [8] 

—I _ I _ I _ I _ I _ 

6  7  8  9  10 

RCpX  10'^ 


20 


Comparison  of  published  heat  transfer  data  for  wire  colls 
and  twisted  tapes. 
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that  pitches  larger  than  about  8.66  had  very  little  influence  on  the  heat 

transfer.  Data  were  reported  as  curves  of  heat  transfer  coefficient  as 

a  function  of  mass  velocity.  The  pressure  drop  data  are  also  given  as 

a  function  of  mass  velocity.  These  pressure  drop  data  were  not  corrected 

for  kinetic  energy  losses.  To  compare  Royd's  results  in  the  form  of 
0 . 4 

NUp/Pr  ■  as  a  function  of  the  bulk  temperature  was  assumed  to  be 

550°F.  These  data  are  given  in  Fig.  4.  The  pressure  drop  data  are  not 
included  in  Fig.  5  as  it  was  not  possible  to  determine  the  necessary 
corrections  for  kinetic  energy  losses. 

Colburn  and  King  [5]  also  investigated  enhancement  effects  of 
twisted  tapes  using  the  same  apparatus  outlined  in  the  preceding  sec¬ 
tion.  Two  full  length,  apparently  snug  fitting  tapes  were  tested; 
however,  data  were  tabulated  for  only  one  of  these  tapes.  Data  were 
taken  for  inlet  air  temperatures  ranging  from  108°C  to  392°C  and  inlet 
tube  wall  temperatures  from  9 . 3°C  to  12.3°C.  As  with  wire  coil  enhance¬ 
ment,  data  are  presented  both  graphically  and  tabularly  as  h  /c  and 

a  p 

corrected  pressure  drop  as  functions  of  mass  velocity  through  the  test 

0  4 

section.  Data,  recalculated  as  Nu^/((Pr)  ’  and  f  as  functions  of 
RCj^  are  included  in  Figs.  4  and  5  for  a  3  ft  long  steel  tape  with  a 
pitch  of  2.67  in  the  2.625  in.  i.d.  tube. 

Studies  on  twisted  tape  enhancement  were  also  conducted  by  Evans 
and  Sarjant  [7]  using  electrically-heated  air  in  an  apparatus  somewhat 
similar  to  that  of  Colburn  and  King.  Tapes  of  2.5  in.  width,  with 
pitches  of  2.33,  3.00,  4.00,  and  4. 67,  were  tested  in  an  8  ft  long. 


3  in.  i.d.  test  section.  The  outside  of  the  tube  was  water  cooled  by 
helically-wound  copper  tubing  soldered  on  the  surface.  Provision  was 


made  for  making  radial  gas  temperature  and  velocity  measurements  at 
several  (unspecified)  axial  locations  so  that  bulk  temperature  and  mean 
velocity  could  be  calculated.  Energy  balances  between  cooling  water 
energy  gain  and  air  energy  loss  were  used  to  validate  data  for  each  run. 
Inlet  air  temperatures  were  approximately  700°F,  while  inlet  Reynolds 
numbers  varied  from  about  3000  to  13,000. 

Results  are  presented  as  curves  of  convective  heat  transfer  coeffi¬ 
cient  and  pressure  drop  corrected  for  kinetic  energy  losses  as 

functions  of  Reynolds  number  for  different  twisted  tapes.  Data  for 

0  4 

twisted  tape  pitches  of  2.33  and  4.67  were  converted  to  Nup/(Pr)  '  as 
a  function  of  Re^  and  are  shown  in  Fig.  4.  Friction  factor  data  for  the 
same  tapes  are  shown  in  Fig.  5. 

Evans  and  Sarjant  attempted  to  separate  the  effects  of  radiation 
and  convection  in  twisted  tape  enhancement.  While  noting  that  their 
analysis  is  not  conclusive  in  this  respect,  they  estimated  the  maximum 
contribution  of  radiation  as  25  percent  of  the  measured  heat  transfer 
coefficient  when  gas  temperatures  above  about  1000°F  were  encountered. 
The  contribution  of  radiation  drops  rapidly  as  gas  temperature  is 
reduced . 

More  recent  tests  of  the  effects  of  twisted  tapes  on  cooled  air 
flows  were  reported  by  Thorsen  and  Landis  [8].  Air  was  supplied  at 
temperatures  up  to  475°F  to  a  1  in.  i.d.  water  cooled  tube.  Twisted 
tape  pitches  were  3.15,  5.17,  and  8.00,  and  the  tube-tape  gap  was  about 
0.003  in.  Heat  transfer  and  friction  factor  data  for  the  smallest  pitch 
tape  are  given  in  Figs.  4  and  5,  respectively.  The  following  empirical 
correlations  were  given  to  describe  the  data. 


Nup  =  0.023[fn^(D,H)lRep“  Pr‘'"  j 


(1  -  0.25 


^  I 

Rt> 


f 


0.046 


Re 


0.2 


12) 


where  fn^CD.H)  and  fn^CDjH)  are  given  in  (8j.  A  major  concern  of  this 

work  was  assessment  of  the  effects  of  variable  properties.  The  tempera 

ture  ratio  correction  factor  was  determined  primarily  from  heating  data 

due  to  the  limited  values  of  (T  /T,  )  which  could  be  attained  in  cooling 

s  b 

The  maximum  value  of  ^  ^  cooling  was  1.05;  however,  the 

scatter  of  the  data  was  at  least  5%. 

Mechanism  of  Lnharuement 

Wire  coils  are  a  popular  turbiilator  insert.  In  the  usual  case 
where  the  coils  are  uniform  in  diameter  and  in  contact  with  the  tube 
wall,  the  wire  acts  .as  a  roughness  element  Ihe  flow  separation  and 
reattachment  are  similar  to  that  pn.dmetl  liv  integral  roughness  ele¬ 
ments  created  by  rolling  or  swaging  The  wires  .are  usually  in  poor 
thermal  contact  with  the  tube,  tiowever.  the  I i n  effect  of  any  of  these 
roughness  elements  is  rather  small.  At  lisu.ii  velocities,  the  flow  is 
predominantly  axial  rather  than  a  swirl  flow  which  follows  thi'  i o i 1 . 

A  heat  transfer  -  momentum  tr.insfei  an.ilogv  s"lution  is  .iv.nlabie  for 
tieliial  repeated  rihs  i |  ,  hut  it  has  iinl  been  tested  against  coil 
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variable  properties  in  cooled  gas  flows  (10],  there  are  no  data  for 
cooling  in  tubes  which  have  repeated  rib  roughness. 

Wire  coils  with  varying  radii  function  more  as  displaced  enhance¬ 
ment  devices,  as  all  but  one  coil  in  each  element  disturb  the  core 
flow  rather  than  the  boundary  layer  at  the  wall.  The  eddies  produced 
in  the  core  are  expected  to  have  some  influence  on  transport  of  momentum 
and  energy  at  the  wall;  however,  the  tf tect  is  probably  less  than  if 
all  coils  were  at  the  wall.  No  models  have  been  proposed  for  this 
situation. 

Twisted  tapes  enhance  heat  transfer  by  creating  a  higher  velocity 
due  to  the  longer  path  length,  by  introducing  a  secondary  flow  due  to 
the  centrifugal  force  acting  on  a  variable  density  fluid,  and  by 
utilizing  the  fin  effect  of  the  tape.  In  the  tests  reported  above,  the 
centrifugal  force  actually  reduces  the  heat  transfer  from  the  air  flow. 
Fin  effects  are  generally  small  because  the  tape-tube  gaps  are  usually 
rather  large  so  that  the  tapes  can  be  easily  inserted  and  removed.  No 
analyses  have  considered  all  of  these  factors;  hence,  empirical  corre¬ 
lations  are  utilized,  e.g.  (8],  which  also  includes  a  correction  for 
temperature-dependent  properties . 


APPARATUS 

An  apparatus  has  been  developed  to  determine  the  thermal-hydraulic 


performance  of  various  turbulator  inserts  with  cooling  of  air.  The 
test  rig  consists  basically  of  a  blower  suitable  to  produce  Reynolds 


16 


riumbers  of  order  10,000,  a  heater  box  to  heat  the  air  to  over  300°F,  and 
the  actual  test  section  where  the  performance  of  the  inserts  is  to  be 
evaluated.  The  cooling  water  jacket  of  the  test  section  simulates 
fire-tube  boiler  conditions  while  serving  as  a  calorimeter. 

Flow  System 

The  U-shaped  flow  system  to  provide  hot  air  to  the  test  section 
is  depicted  in  Fig.  6.  Airflow  is  provided  by  a  centrifugal  blower  that 
delivers  600  SCFM  at  a  head  of  6  in.  of  water.  After' leaving  the  fan, 
the  air  passes  through  flow  straighteners  and  a  flow  development  length 
of  10  ft  before  entering  an  orifice  meter.  Eight  feet  downstream  from 
the  orifice  meter,  the  air  passes  through  an  elbow  before  entering  the 
heater  box. 

The  heater  box  was  fabricated  from  22  gage  sheet  steel  in  three 
parts:  an  inlet  rectangular  diffuser  section  varying  in  cross  section 

from  6  in.  x  6  in.  at  the  entrance  to  21  in.  x  il  ]/4  in.  over  a  length 
of  72  in.  The  flow  is  at  low  velocity  as  it  enters  the  8  in.  long  sec¬ 
tion  containing  a  Chromalox  Process  Air  Duct  Heater,  Model  ADHT-OIOXX, 
rated  at  10  KW,  240  V,  3  phase.  The  box  cross  section  is  reduced  from 
21  in.  X  11  1/4  in.  back  to  6  in.  x  6  in.  over  a  length  of  50  in.  down¬ 
stream  of  the  heater.  Each  part  of  the  box  was  fabricated  in  two  halves 
and  joined  by  welding  along  the  center. 

Three  autotransformers  are  connected  in  a  208  V  Delta  circuit  to 
provide  the  electrical  input  to  the  heater.  High  temperature  insulating 
gasket  material  is  used  to  prevent  direct  contact  between  the  heater  and 


the  duct.  An  over-temperature  protector  is  included  in  the  electrical 
circuit  to  cut  off  electrical  power  to  the  heater  once  a  predetermined 
temperature  is  reached.  A  motorized  mixer  is  provided  to  minimize 
thermal  stratification  within  the  heater  duct. 

The  hot  air  coming  out  of  the  heater  box  passes  through  a  bend 
into  a  developing  length  of  20  ft  before  entering  the  test  section. 
After  passing  through  the  test  section,  the  air  is  exhausted  into  a 
5  ft  long  pipe  section  before  discharge  outside  the  building.  The 
heater  box  and  piping  leading  into  and  out  of  the  test  section  are 
insulated  with  calcium  silicate  blocks  capable  of  withstanding  tempera¬ 
tures  of  more  than  1500°F.  These  blocks  are  covered  with  an  additional 
layer  of  glass  fiber  to  reduce  the  heat  loss  to  a  minimum. 

The  test  section  is  insulated  simply  with  glass  fiber  since  the 
outside  temperature  of  the  calorimeter  is  less  than  150°F.  Precautions 
were  also  taken  to  insulate  the  hypodermic  tubing  carrying  thermocouple 
leads  for  measuring  temperatures. 

Test  Section 

The  test  section  shown  in  Fig.  7  is  based  on  a  3  in.  o.d.  carbon, 
drawn-on-mandrel  steel  tube  with  i.d.  of  2.675  in.  to  provide  a  snug 
fit  with  the  inserts  provided.  The  tube  is  71  3/4  in.  long  and  is 
flanged  at  the  ends  to  facilitate  installation  and  removal  of  inserts. 
Each  end  of  the  tube  is  fitted  with  three  pressure  taps  spaced  equally 
around  the  circumference.  Tube  wall  temperatures  are  measured  with 
chrome  1 /a  1 ume 1  thermocouples  fitted  into  axial  slots  milled  part  way 
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through  the  tube  wall.  Five  axial  locations  for  these  thermocouples 
were  chosen  to  correspond  with  the  water  calorimeter  cooling  sections 
described  below.  Three  thermocouples  were  installed  at  each  axial 
location,  spaced  equally  around  the  tube  circumference.  The  thermo¬ 
couples  were  made  of  welded  junctions  and  were  inserted  into  the  slots 
through  small  tubes  passing  under  the  calorimeter  cooling  sections. 

The  water  calorimeter  for  cooling  the  gas  and  conducting  energy 
balances  is  made  of  1/4  in.  i.d.  copper  tubing  wound  around  the  tube. 
Prior  to  installation,  the  outside  of  the  steel  tube  was  cleaned  and 
tinned.  Then  the  copper  tube  was  soldered  to  the  outside  with  a  low 
temperature  (430°F)  silver  bearing  solder.  All  State  No.  430  with  Duzall 
liquid  flux.  Heating  was  by  a  torch;  to  facilitate  rotation,  the  tube 
was  installed  in  a  lathe  and  the  copper  tubing  was  soldered  to  the 
steel  tube  as  it  was  being  tightly  wound. 

The  calorimeter  tubing  was  connected  as  four  separate  subsections 
in  series,  with  water  temperature  measured  at  the  inlet  and  outlet  of 
each  subsection.  This  feature,  together  with  the  wall  thermocouples 
installed  as  described  above,  permitted  the  determination  of  four  sec¬ 
tional-average  gas  side  heat  transfer  coefficients. 


Instrumentation 

The  basic  measurements  needed  are  the  air  and  water  mass  flow  rates 
and  the  various  air  and  water  bulk  temperatures.  Also  required  for 
calculation  of  the  heat  transfer  coefficients  are  the  wall  temperatures. 


Additionally,  the  pressure  drop  is  required  to  calculate  the  friction 
factor . 


The  air  flow  is  controlled  by  a  blast  gate  at  the  blower  outlet  and 
passes  through  flow  straighteners  before  entering  the  section  upstream 
of  the  ASME  orifice.  The  pressure  drop  across  the  orifice  is  monitored 
by  a  differential  pressure  transducer,  which  is  connected  via  a  Scanivalve 
to  a  piezometer  ring  attached  to  the  multiple  pressure  taps.  The  same 
system  was  used  to  measure  static  pressure.  The  mass  flow  rate  was 
found  using  standard  ASME  procedures  and  equations.  Also  measured  were 
the  static  pressure  at  the  test  section  inlet  and  pressure  drop  across 
the  test  section. 

The  temperatures  measured  are  the  wall  temperatures  along  the 
test  section,  inlet  and  outlet  air  temperatures,  water  temperatures  at 
the  inlet  and  outlet  of  each  segment  of  calorimeter  tubing,  and  the 
room  temperature.  Thermocouples  are  made  of  28  ga  chromel/alumel  wire, 
referenced  to  a  150°F  Whitaker  reference  junction.  Both  temperature 
and  pressure  signals  were  recorded  digitally  by  the  Heat  Transfer 
Laboratory  Data  Acquisition  System  [11].  Finally,  the  water  flow  rate 
was  measured  directly  using  a  container  on  a  sensitive  balance  and  a 
stopwatch . 

Some  of  the  problems  anticipated  in  the  apparatus  were  air  leaks, 
stratification  of  air  flow  in  the  heater  box,  axial  conduction  of  heat 
to  the  test  section,  and  radiation  errors  in  temperature  measurements. 

The  absence  of  air  leaks  was  confirmed  by  comparing  air  flow  measured 


by  the  orifice  downstream  of  the  blower  with  flow  measured  at  the  exhaust 


22 


(past  the  test  section)  using  a  pitot  tube  and  taking  a  velocity  profile. 
The  difference  was  less  than  one  percent. 

As  noted  earlier,  a  mixer  was  introduced  into  the  heater  to  heat 
the  air  more  uniformly.  This  avoided  early  tripping  of  the  power  supply 
by  the  over-temperature  protector,  whose  sensing  thermocouple  is 
located  in  the  middle  of  a  heater  element  at  the  downstream  last  row  of 
vertically  heated  elements.  To  minimize  the  axial  conduction  from  the 
test  section,  two  teflon  spacers,  2  1/2  in.  and  1  1/2  in.  thick,  were 
placed  at  the  inlet  and  outlet  of  the  test  section,  respectively.  Four 
thermocouples  were  inserted  on  either  side  to  measure  the  temperature 
drop  across  the  spacers  so  that  the  axial  conduction  could  be  estimated. 
Regarding  the  influence  of  radiation  on  temperature  measurements,  cal¬ 
culations  showed  that  the  errors  were  very  small  and  well  within  the 
experimental  error  margin. 


Procedure 


Testing  and  Data  Acquisition 

After  starting  the  blower  and  the  water  pump,  the  electrical 
heater  was  switched  on.  The  over-temperature  controller  was  set  pre¬ 
viously.  Initially,  the  air  flow  was  kept  low  so  that  the  system  could 
warm  up  faster.  Normally,  the  time  taken  for  the  air  temperature  at 
the  test  section  inlet  to  reach  300°F  was  approximately  three  hours. 

The  stirrer  raised  the  temperature  of  the  air  considerably  (sometimes 
up  to  40°F)  through  uniform  mixing  around  the  heater.  The  flow  rate  was 
then  adjusted  for  the  run  requirements.  When  steady  state  was  reached, 
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the  temperature  drop  from  the  heater  box  to  the  test  section  was  less 
than  25°F.  Whenever  the  flow  rate  (Reynolds  number)  was  subsequently 
changed,  it  took  about  twenty  minutes  to  attain  steady-state  conditions 
on  both  air  and  water  sides.  The  water  flow  rate  was  set  at  about 
0.2  gpm  which  represented  a  good  compromise  between  an  accurately 
measurable  temperature  rise  and  a  nearly  isothermal  tube  wall. 

Once  the  steady  state  was  reached,  the  reading  of  various  tempera¬ 
tures  and  pressures  could  be  done  in  about  one  minute  through  the  data 
acquisition  system.  The  system  used  was  identical  to  the  one  described 
in  Reference  11,  with  the  exception  of  the  digital  voltmeter,  which  was 
a  more  recent  model.  A  program  was  written  which  enabled  the  system  to 
take  all  the  required  measurements  as  well  as  reduce  them  to  relevant 
information.  The  program  flow  chart  and  listing  are  given  in  Appendix  B. 
Data  Reduction 

The  computer  was  instructed  to  scan  inlet  air  and  outlet  water 
temperatures  for  steady  state.  At  steady  state,  all  thermocouples  and 
previous  transducer  readings  were  read,  quantities  were  calculated,  and 
the  heat  balance  was  checked: 

q  -  q  =  q  (3) 

^w  ^c  a 

Agreement  to  within  5%  was  required  to  validate  a  run.  (Most  energy 
balances  agreed  to  within  2-3%. )  If  acceptable,  the  data  reduction  was 
concluded  and  the  results  were  printed  out. 

Average  heat  transfer  coefficients  were  determined  for  each  calori¬ 
meter  segment  of  the  test  section.  Starting  at  the  first  segment,  the 

energy  added  to  the  water  is  obtained  from  q  -q  =mc  (T„-T,) 

^w  ^c  w  pw  w2  wl 


where  is  the  estimated  heat  transfer  to  the  water  by  conduction  from 

the  upstream  ducting,  q^  =  k^A(AT^/AX).  Assuming  that  ~  " 

me  (T  ,  -  T  T  „  can  be  calculated.  This  involves  iteration  since 
pa  al  a2  a2 

the  air  specific  heat  is  dependent  on  temperature.  The  average  surface 
temperature,  T^ ,  is  taken  as  the  average  of  the  two  average  temperatures 
for  the  two  measuring  stations  in  each  section,  that  is,  the  linear 
average  of  6  readings.  The  heat  transfer  coefficient,  h  ,  is  then 

3 

obtained  from  the  usual  expression 
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The  process  is  repeated  using  the  outlet  air  temperature  at  one  section 
as  the  inlet  temperature  for  the  next  section.  The  second  and  third 
segments  have  no  conduction  correction;  however,  there  is  heat  gain  by 
the  coolant  in  the  fourth  section. 

In  calculating  the  Reynolds  number,  no  allowance  was  made  for  the 
blockage  of  turbulators.  The  average  Fanning  friction  factor  is  given 
by 
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where  the  measured  pressure  drop  is  corrected  for  the  momentum  change. 
Again,  the  cross-sectional  area  is  that  of  the  empty  tube  without  con¬ 


sideration  of  turbulators. 


Uncertainties  in  Nusselt  number,  Prandtl  number,  Reynolds  number 
and  friction  factor  were  estimated  from  a  typical  run  at  Re  =  12,500. 


The  uncertainties  given  below  are  typical  of  the  uncertainties  that 
can  be  expected  in  other  test  runs. 


Term 

Uncertainties 

Nu 

±5% 

Re 

±3% 

Pr 

±3.5% 

f 

±5.3% 

Error  analysis  details  and  sample  calculations  are  given  in  Appendix  D. 


RESULTS  AND  DISCUSSION 
Results 

Before  starting  to  evaluate  the  performance  of  the  inserts,  the 
experimental  apparatus  was  tested  for  its  performance  as  a  plain  tube 
without  using  inserts.  Flows  for  Reynolds  numbers  from  6,000  to  25,000 
were  tested.  Heat  transfer  and  friction  characteristics  were  determined 
for  nominal  inlet  air  temperatures  of  330°F.  Errors  in  the  energy  balance 
were  under  five  percent,  with  most  balances  falling  within  two  or  three 
percent.  The  results  compared  favorably  with  the  well-accepted  equation 
of  Dittus  and  Boelter  [12], 

Nu^  =  0.023(Re  )°-®(Pr)°-^ 


(6) 
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as  shown  in  Fig.  8.  Experimental  data  for  Fig.  8  are  given  in  tabular 
form  in  Appendix  C. 

Test  runs  were  conducted  over  a  range  of  Reynolds  numbers  for  the 
three  inserts.  The  reduced  data  were  plotted  as  curves  of  NUj^  vs. 
and  f  vs.  Re^  as  shown  in  Figs.  9  and  10.  The  Type  A  and  Type  B  inserts 
produce  large  increases  in  heat  transfer  coefficients  over  the  smooth 
tube  values.  Coefficients  are  increased  a  lesser  amount  with  the  Type  C 
insert . 

Experiments  were  performed  for  horizontal  and  vertical  orientations 
of  the  Type  A  and  Type  B  inserts;  the  data  are  compared  for  each  insert 
in  Figs.  11  and  12.  Performance  differences  for  the  two  orientations 
are  small  except  for  the  Type  B  turbulator  at  lew  Reynolds  numbers  where 
a  slight  decrease  is  found  for  the  horizontal  orientation. 

Discussion 

The  Type  A  and  turbulators  produce  about  150%  increase  in  Nusselt 
number  over  a  plain  tube  while  the  Type  B  turbulator  provides  about  12_% 
increase.  The  two  turbulators  are  actually  quite  close  in  performance 
considering  the  measurement  uncertainty.  Also,  operating  conditions  in 
actual  boilers  may  make  any  difference  of  negligible  importance.  The 
Type  C  insert  showed  an  increase  in  Nusselt  number  about  60  percent 
above  that  of  the  plain  tube.  Some  visual  observations,  made  with  a 
plexiglass  tube  and  smoke  generator,  indicated  that  the  mixing  was 
considerably  better  with  the  Type  A  and  Type  B  turbulators  than  with 
the  Type  C  insert. 


NUSSELT  NUMBER,  Nu^-fd/D) 


Fig.  8.  Experimental  data  for  smooth  tube 


Fig.  11.  Comparison  of  heat  transfer  data  for  vertical  and  horizontal 
orientation  of  insert  Type  A. 
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WTifr,  chf  'VTP  inseits  Wfre  positioned  horizontal  iy ,  the 

Nusselt  number  differences  over  the  range  of  the  Reynolds  number  tested 
are  within  experimental  error  except  for  the  Type  B  insert  between 
Reynolds  numbers  of  6,000  and  10,000.  The  difference  was  between  15  and 
8%  up  to  a  Reynolds  number  of  10,000.  At  higher  Reynolds  numbers,  the 
difference  is  negligible.  It  may  be  that  in  the  horizontal  orientation 
at  lower  Reynolds  numbers,  the  mixing  is  not  complete  so  that  the 
hotter  air  travels  along  the  upper  part  of  the  tube  in  a  partially 
stratified  flow.  However,  no  detailed  data  were  taken  to  confirm  this 
speculation . 

A  review  of  the  experimental  data  given  in  Table  11  lor  Reynolds 
numbers  of  10,700  and  25,000  indicates  more  specifically  the  increases 
in  Nusselt  number  and  friction  coefficient  over  plain  tube  data. 

The  friction  coefficient  for  the  Type  A  insert  is  12.1  times  that 
of  a  plain  tube  at  a  Reynolds  number  of  about  10,000;  the  increase  is 
13.8  times  at  a  Reynolds  number  of  25,000.  The  Type  B  insert  has  f 
increases  of  11.0  and  12.4  times  for  the  same  Reynolds  numbers,  while 
the  Type  C  insert  has  increases  of  2.7  and  2.4  times,  respectively, 
for  these  Reynolds  numbers.  The  friction  and  heat  transfer  character¬ 
istics  of  the  Type  C  insert  are  both  '■-•mparable  to  data  noted  in  the 
Literature  Review. 

When  evaluating  the  relative  merits  of  pressure  loss  cost  tor  a 
particular  heat  transfer  coefficient  i  . crease,  it  should  be  borne  in 
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Table  II.  Nusselt  number  and  friction  coefficient  increases 


Insert 

Rep 

Plain 

Plain 

Nuj^ 

f 

% 

f 

Type 

% 

f 

Ratio 

Ratio 

Type  A 

10,700 

39 

.00785 

100.4 

.095 

2.57 

12.1 

25,000 

76.5 

.00630 

180 

.087 

2.35 

13.8 

Type  B 

10,700 

39 

.00785 

87.5 

.086 

2.25 

11.0 

25,000 

76.5 

.00630 

160 

.078 

2.09 

12.4 

Type  C 

10,700 

39 

.00785 

62.5 

.021 

1.60 

2.7 

25,000 

76.5 

.00630 

110.1 

.0151 

1.44 

2.4 

mind  that  heat  transfer  can  be  enhanced  considerably  by  merely  increas¬ 
ing  the  air  velocity  with  a  concomitant  increase  in  pressure  loss. 

This  situation,  noted  many  years  ago  by  Colburn  and  King  15],  is  a 
possible  alternative  to  augmentation  by  use  of  inserts.  Said  differ¬ 
ently,  if  a  given  increase  in  pressure  loss  is  acceptable,  it  is 
possible  to  obtain  the  same  augmentation  by  increasing  the  velocity  (and 
the  Reynolds  number)  to  obtain  a  higher  convective  coefficient.  However, 
this  option  is  generally  not  available  in  a  fire  tube  boiler. 

Finally,  it  should  be  noted  that  a  "fine  scale"  variation  of  Nusselt 
number  as  a  function  of  Reynolds  number  may  be  discerned  for  each  sub¬ 
section  of  the  test  section,  even  though  all  data  fall  within  the  con¬ 
fidence  bounds  given.  Sufficiently  detailed  tests  to  investigate  and 
comment  on  this  second-order  phenomenon  have  not  been  included  in  this 
study.  It  was  determined,  however,  that  the  temperature-dependent 
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property  correction  [Eqs.  (1)  and  (2)]  did  not  significantly  reduce  the 
scatter . 


CONCLUSIONS 

The  results  show  that  connraercial  turbulator  inserts  provide  con¬ 
siderable  enhancement  of  gas-side  convective  coefficients  in  fully- 
developed  tube  flow.  It  may  be  possible  to  reduce  fuel  consumption 
through  the  use  of  inserts,  but  improvement  in  gas-side  convection  co¬ 
efficients  does  not  result  in  ?  one-to-one  ratio  improvement  in  boiler 
efficiency.  Rather,  augmentation  of  gas  side  coefficients  is  only  one 
factor  among  many  items  such  as  fuel/air  ratio  adjustment,  tube  cleanli¬ 
ness,  and  air/gas  system  flow  resistance  that  help  reach  goals  of 
reduced  fuel  consumption  and  improved  boiler  efficiency.  It  should 
further  be  noted  that  calculation  of  boiler  efficiency  with  inserts 
should  consider  changes  in  the  water  side  flow  distribution  (and  heat 
transfer  coefficient)  brought  about  by  the  large  increase  in  gas-side 


heat  transfer  coefficient. 
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Brock  "Fuel  Saver"  Turbulator 
Fuel  Efficiency,  Inc. 

131  Stuart  Avenue 
P.O.  Box  48 

Newark,  New  York  14513 
(Manufacturer) 

Smick  Turbulator 
American  Fuel  Economy,  Inc. 

Rt.  1  Canal  Rd.  E. 

Ottawa,  Illinois  61350 

(Manufacturer.  Direct  sales  to  boiler  manufacturers) 

AEI  Fire  Tube  Turbulator 
Allied  Energy  International 
8601  Wilshire  Blvd.,  Suite  1002 
Beverly  Hills,  California  90211 
(Retrofit  distributor) 

Powermaster  Turbulator 
Powermaster  Corp. 

P.O.  Box  231 

Middletown,  Pennsylvania  17057 
(Manufacturer) 


A 


A5 


U  :  dim  I  J  5)  ,PS  1  I'll  ,  x$l  ]  5|  ,  X(  17)  ,D|471  ,B  (  11)  ,FI  11]  ,GH1  I 

1:  dim  C147|,R(7|,E|7] 

2:  wtb  6, 32  ,  32 , 32 ,  32  ,  32  ,  32  ,  32,  32,32 
3  :  wtb  6,27,77 

4:  dsp  "Set  printer  at  t<3p  of  pa9e“;sto 
3;  wtb  6,27,84 

6:  wtb  6 , 27 , 7 0  ,  1  n t (  1 0 56/6 4)  , i nt ( 10 56) 

7:  wtb  6,27,76, int(960/64) ,int(960) 

8:  ent  "Operator?",AS 
9 :  ent  " Da  te  ? " , B  S 
10:  ent  "Type  of  Pun?",XS 
11:  fmt  l,"C",fz2.0,"E",2 
12:  fmt  2 , "F1R3HOA1T1M3" 

13:  0-C;wrt  7d9.1,C;wrt  722,2 
14 :  for  1=1  to  3 

15:  beep;dsp  "Set  scannivalve  on  clianne  1 "  ,  I ;  s  tp 

16:  if  I  =  3:ent  "Max.  Pres,  reached?  Yes= 1 : No=0"  ,N ; i f  N=0;qto  18 
17:  if  1=3, -ent  "Pres,  in  i n . H20? "  ,D [ 3 2 1 1 ] ;9 to  20 
18:  0-K;for  N  =  1  to  10;rcd  722 , G ;G  +  K-K ; nex t  N ; K/1 0-C I  3 2 U I 
19:  -1.13  9C(33)+1.13  9C|32t-I|*D132H| 

20:  next  I 

21:  beep;ent  "Atmospheric  pressure , in . Hq " ,C [  39)  ;C ( 39 | »D I  3 9 ) 

22:  fmt  3 , " FI R2A IHIT 1M3 " 

2  3:  32*C;wrt  709.1,C,-wrt  72  2.3 

24:  0-K;for  N=1  to  10;red  722 , G;G+K-K ; nex t  N ; lOOOK / 1 0-C ( 3 2 ) 

25:  14^  1  195+4  4  .  1465C1  32) +  .05480  32)  *2-. 00510  32)  ‘3*D|  321 
26:  (1.3259D|39|-.097527D(34|)/(D(32|+459.67) -X) 1] 

27:  2.6866e-3(Dl 32) +459.67) '1.5/(D[32) +658. 37) -XI 2) 

28 :  15000-RI 1) ; 25000-Rl 2) ;35000-R(3) ; 50000-RI 4) ;75000-R[ 5| ; 1 000 00-R [ 6 ) 

29:  150000-R(7i  ; .026  12-E (1  I ; .02601-E(2) ; .02596-E(3) ; .  02591-E (4 ) ; .02585-E (51 
30:  .02582-E!61  ;  .02579*E|7) 

31:  E(11*B;P(1|*R 

32:  12932.9505*B*V(XI1I*D( 351)-M 

33:  12. 2525M/X( 21 -S 

34:  if  abs  ( S-R)  >.  1 0 ;  q  to  40 

35:  for  1=1  to  6 

36;  if  S>R( 1) ;gto  38 

3  7:  no  X  t  I 

38:  E[Il-(3-F(Il)*(E(I|-F(I+l))/(R(I+l)-R(I))-B 
39:  S-R;gto  32 
40:  M-X131 

41:  dsp  "Mass  flow,lbm/hr=  ",X(3);stp 

42:  ent  "Is  Mas:,  flow  rate  OK?  Yes  =  l ;  No=0"  ,N 

4  3:  if  N=0 ;g  to  11 

44:  dsp  "Set  Air  Tempe r a t ure "  ;s tp 

45:  fmt  3  ,  "  F1F2A1.H  1T1M3  "  ;  f  mt  1 ,  "C"  ,  f  z2 . 0  ,  "F"  ,  z  ;  wr  t  722.3 
46:  24-C;wtt  709. 1,C 

47;  red  72 2 ,G ; lOOOG-G; 149 , 1 195+4 4 . 146 5G+ . 0548G ' 2- . 0 0 5 IG ‘ 3-G 
48:  d.sp  "Outlet  Water  Temp=  ",G;wait  1000 
49:  d.sp  "Inlet  Air  T8mo=  ",F;wait  1000 
50:  1-C;wrt  709. 1,C 

51:  red  722  , F ; lOOOF- F ;  14 9 . 1  1 95  +  44 . 14 6 5F+ . 0 548 F" 2- . 00 5 1 F ' 3 - F 
52:  dsp  "To  Record  Da ta ; STOP, cont 57 , EXC" ; wa i t  1000;gto  46 
53:  fmt  3 , " F 1 R2A  IHlT 1 M 3 " ; f mt  1 , "C” , f z2 . 0 , "E " , z ; wr t  722.3 
54:  for  1=1  to  32;l--:wrt  709. 1,C 

55:  0-K;for  N=1  to  10;red  722 , G;G+ K-K ; ne x t  N ; 1 OOOK/1 0-C ( I ) 

56  :  149. 19950354  +  44. 146471  3C[I)  +  . 0547791 6C (I)' 2-. 00 50 9771C[I)'3-D1I] 

5  7:  no  X  t  I 

58:  for  1=1  to  4;I+32-C;wrt  709. 1,C 

59:  0-K;for  N=1  to  10;red  722 , G;G  +  K-K ; nex t  N ; 1 OOOK/ 1  0-C (  I +4 3 ] 


-1 


•J 


A 


1 


.1 


Fig.  B  -  2.  Data  acquisition  program 


H9.  199'i*-44.  1465CI  I  +4  3J  +  .  0  548C1 1  +4  3  |  “  .  0  0  5 1C  I  I  4  :i  J  '  3  -D  |  I -►  4  3  J 

ne  X  t  I 

£iPt  2  ,  "F1R3HOA1T1M3" 

0'»C;wi't  709.1,C;wrt  722.2 
foe  1=1  to  6 

beep;dsp  "Set  Scannivalve  at  Channel  ",I;sto 

if  I  =  3;ent  "Max.  Pres,  reached?  Yes= 1 : No  =  0 "  ,N ;  i  f  N=0;qto  68 

if  I  =  3;ent  "Pres,  in  i n . H20? "  ,D  I  32 1 1 1  ;q to  70 

0*K;£or  N  =  1  to  10:red  72  2  ,  G  iG  +  K-'K  ;  ne  x  t  N  ;  K  /  1 0*C  1  1  +  3  2  I 

-1.139Cf33i+1.139CiI+32)-DII+32| 

no  X  t  I 

ent  "Pnter  Atmos. Pres,  in .Hq",C(39);C139|*D139] 

ent  "Water  flow  ratein  Ibm/hr ? " , D ( 4 2 | 

n [ 42)  *0  1 40 i -D( 4  11  *D( 4  31 

(D (24 1 +D[ 251  eD [261  +D( 27|  ) /4  *X(4  I 

(D|281+D(291+D[30|+D|31|)/4*X(51 

X(  41  -Xi  5)  »X| 61 

(X  1  4  I  +X  1  51  )  /2-*r0 

.2231+3.42e-5(r0+459.67)-2.93c-9(r0+459.67)'2*X|71 

(1.3259Di39|-.097327D|34|)/(D132|+459.67).Xil! 

2. 6866e- 3(D|  321 +4  59 . 67) "1 . 5/(Dl 32  ) +658 .  37) ♦X i 21 
f.  I  1  I  :  P  1  1  1  •  P 

12932. 9505 *B+((Xil)*D[35|)*M 

12, 2525M/XI 2| -S 

if  abs ( S-P ) < i 0 ; q to  9J 

for  1  =  1  t  o  6 

if  S  '  P  I  I  1  ;  1  .Ti  p  2 

next  I 

r|I)-(S-R(Il  )*(FlIi-f:|I  +  l|)/(P|l+l|-RlI  l)-B 
S  ■*  P ;  q  t  o  8  2 

M*Xl  31  ;B-X(  1  3|  ;S-X|  151 
XI  31*X|71 *X|61-X|8| 

DI4  31  (Dl  17)  -D|  16)  )  -Xi  10) 

D1  421  (D(  191  -D|  18  1  )  -X  r.  2) 

D( 4 1  I  (D( 211  -D( 20  I ) *X|  14  1 
0(401  (Dl  231 -D( 221  ) -Xl  16) 

X(  101  +X1  12]  +X( 14 1+Xl 16 1-XI 17) 

.759255(0147] -D(46| )*Xf91 
1.071844(D(44|-D[45i )*X(11| 

XI  41  *r9:Xl  3|  TlOjXl  161  -X(9|*rn;  .971429T  16 
(Dl  1|  +D|  21  +D1  3|  )/'3*A  I  3|  -rl4  ;  (D|4  )+D(  5|  +D|  6)  )/3-A  (4  1  -.rl5 
cll  'nirm'(r9,rl0,rll,rl2,rl3,rl4,rl5,rl6,rl7,rl8,rl9,r20,r21,r22) 
r  11-A  I  1  I ; r 12-AI 2 !  ; r 16*A1 51  ; r 17*A (61  ; rl 8-A( 7] ; r 19*A (8) 
r20-A(91 ; r21-Al  10  !  ;r 22  -Al 1 1  I ;A( 1 ) ♦r9;Xl  141  *r 1 3;  1  .  0  30  30  3-rl6 
AI4)*B13)-rl4;(D171+D|8)+D(91)/3-B(4)*rl5 

cll  'ntrm'(r9,rl0,rll,rl2,rl3,rl4,rl5,tl6,rl7,rl8,rl9,r20,r21,r22) 
r 11-B (1  I ; r 12-B (2|  ; r 16-B (5)  ; r 17-B  (6|  ; rie.B ( 71 ; r 19-B (8  I 
r20-B  (9|  ;  r-?  1  -  E  (  10  1  ;  r  22-B  (  1 1  1  ;B  I  1  1  ♦r9;Xl  12|  *r  1  3;  1  .007407-  1-16 
B141-F(3)-rl4;(DllUi+D|ll)+D|12))/3-F(4)-rl5 

cll  'nirm'(r9,rl0,rll,rl2,[13,rl4,rl5,rl6,rl7,rl8,tl9,r20,r21,ri2) 
rll-F[lJ;rl2-F|2|;rl6-F|5|;rl7-F|61:rl8-Fl71;rl9-F181 
r20-F[91;r21-F(101;r22.F|lll;F|l)-r9;X(10|-X|lll-rli;.992701-rl6 
F (4 1 -G  1  31 -r 14 ;  (D I  1  1  I +D|  14 )+D| 151 ) /3-GI 4  I -r  15 

c-11  nirm'(r9,rl(!,rll,rl2,rl3,tl4,rl5,rl6,rl7,rl8,rl9,r20,r21,r22) 

r  1  1  *G  1  1  :  ; r  U  C i 2 i ; r 16-GI 5 1 ; t 1 7-G(6I  : r 18-G (7)  ; r 19-G  |8 I 

r  20  -G  :  9  i  ;  r  -  G  f  i  n  i  ;  r  2  2  -G  (  1  1  1 

d  j)'  "Fr.-‘rqy  Balance  =  "  ,  X  1  8  1 /X  |  1  7  ]  ;  s  t  p 

wr  r  6,"Dat>  :  ",liS,"  "."Operator:  ",AS 

wi r  6, "Type  of  Pun:  ",XS 

w:t  6," -  -  - - - - - - - - - - - 


B  -  2.  Continued. 
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;  20 

wr  t 

6  ^ "  Raw  Data" 

„ 

1  2  I 

wr  t 

122 

wr  t 

6  ,  "  Temperature 

Data" 

12  3 

wr  t 

6, "Local  ion  Scanner 

Read  rnq 

Te  moe  ra  t  ure ' 

1  2  4 

wr  t 

6,"  Channel 

Vc)  i  t  s 

deg . F" 

125 

wr  t 

126 

f  mt 

5,2fl2.6 

127 

wr  t 

6. 5, "Test  Sec. wall  l-tOp,C(l) 

1 

",C|  1)  ," 

",  D I  1  1 

126 

w  r  t 

6. 5, "Test  Sec. wall  l-ns,Clll 

2 

",C|21  ," 

", D 1 2  1 

129 

wr  t 

6. 5, "Test  Sec. wall  l-fs,C(3I 

3 

",CI 31  ," 

"  ,  D  1  3  i 

130 

wr  t 

6. 5, "Test  Sec. wall  2-top, C|4) 

4 

",C|4  1  ," 

",  D 1 4  1 

131 

wr  t 

6. 5, "Test  Sec. wall  2-ns,C(5) 

5 

",CI 51  ,  " 

",  D  1  5  1 

132 

wr  t 

6. 5, "Test  Sec. wall  2-fs,C(6| 

6 

",C!61 ," 

", D 1 6 1 

133 

wr  t 

6. 5, "Test  Sec. wall  3-top,C|7| 

7 

",C[7| 

",  D 1 7 1 

134 

wr  t 

6. 5, "Test  Sec. wall  3-ns,C18| 

8 

",C(8i  ," 

"  ,  D  [  8  1 

135 

wr  t 

6. 5, "lest  Sec. wall  3-fs,C19| 

9 

" ,C 1 9 1 , " 

"  ,  D  (  9  1 

136 

w  r  t 

6. 5, "Test  Sec. wall  4-top, C|101  10 

",C1 101 , " 

"  ,  D  1  1  0  i 

137 

wr  t 

6. 5, "Test  Sec. wall  4-ns,Clll) 

11 

"  ,  C  I  1  1  1  ,  " 

", D  1 1 1 1 

138 

wr  t 

6. 5, "Test  Sec. wall  4-fs,C|12| 

12 

",C| 12]  ,  " 

",0112) 

1  3  9 

wr  t 

6. 5, "Test  Sec. wall  5-top, 0(13)  13 

",C|  13)  ,  " 

", D  1  1 3 1 

140 

wr  t 

6. 5, "Test  Sec. wall  5-ns,Cil4) 

14 

",C(  14  1  ,  " 

",D|141 

14  1 

wr  t 

6. 5, "Test  Sec. wall  5-fs,C115i 

15 

",C|15i  ," 

",D|151 

142 

wr  t 

6. 5, "Water  I nl e t-4 ,C ( 16 1 

16 

",C(16)  ," 

", D1 16] 

143 

wr  t 

6. 5, "Water  Out le t-4 ,C | 1 7| 

17 

",Cil7i  ," 

",D|17) 

144 

wr  t 

6. 5, "Water  I nle t- 3 , C | 1 8 1 

18 

",CI 181  ,  " 

",  D 1 18] 

145 

wr  t 

6. 5, "Water  Out le t-3 ,C( 19  I 

19 

",C( 19) ,  " 

",  D  I  1 9 1 

146 

wr  t 

6. 5, "Water  I n le t- 2, C ( 20 ) 

20 

",CI20)  ,  " 

",D|201 

147 

wr  t 

6. 5,  "Water  Outlet-2, C(211 

21 

",C(21|  ," 

"  ,  D  [  2  1  ] 

148 

wr  t 

6. 5,  "Water  Inlet-1, C|22) 

22 

",C(221 ," 

",D122) 

149 

wr  t 

6. 5, "Water  Out le t- 1 , C ( 2 31 

23 

",C|231 ," 

",D(23I 

150 

wet 

6. 5,  "Alt  In let-1, 0  124) 

24 

",C124) ," 

" ,D|24) 

1  51 

wr  t 

6. 5,  "Air  Inlet-2, Cf25) 

25 

",C(25| ," 

",D|251 

152 

wr  t 

6 . 5 , " A 1  r  I  nle  t-3 ,C ( 26  1 

26 

",C!26) ," 

",D[26) 

153 

wr  t 

6. 5, "Air  Inlet-4, C(27| 

27 

",C|271 ," 

",n|27I 

154 

wr  t 

6. 5, "Air  Ou t le t- 1 , C ( 28 1 

28 

",C(281  ,  " 

",D|281 

155 

wr  t 

6. 5, "Air  Oi;  t  le  t-2  ,C  (  29  ) 

29 

",C[29| ,  " 

",D129) 

156 

wr  t 

6. 5, "Air  Outlet-3, C[30) 

30 

",C(30|  ," 

",D|30) 

157 

wr  t 

6. 5, "Air  Outlet-4, C131I 

31 

",C|3]|," 

", D 1 3 1 1 

158 

wr  t 

6. 5, "Room  Temperature, Cl 32j 

32 

",C|32) , " 

",  D 1 3  2 1 

159 

wr  t 

6,"  Pressure 

Data" 

160 

wr  t 

6  ,  "  - 

161 

wr  t 

6, "Location  Scannivalve  Readinq 

Pressure" 

16  2 

wr  t 

6 , "  Channel 

Volts 

1  n  .  H  2  0  " 

16  . 

wr  t 

^  ,  - - - - - 

164 

wr  t 

6  .  5  ,  "Atmosphere  ref.,C|33| 

1 

"  ,CI 3  3]  ,  " 

"  ,  D  (  3  3  1 

!6' 

w  r  t 

6  .  5  ,  "St  .  inlet  h . h. , C( 34 ) 

2 

",C( 34 1 , " 

",  0  1  3  4  1 

165 

wr  t 

6. 5, "dp.  orificc,C(35i 

3 

",C(  351 , " 

",D(351 

167 

wr  t 

6. 5, "St.  ds . or  1 f 1 cc ,C 1 36 1 

4 

",C! 361 , " 

",  D  I  3 6  1 

168 

wr  t 

6. 5, "dp.  test  sec 1 1 on , C ( 3 7 1 

5 

",C| 371  ,  " 

"  ,  D  1  3  7  ! 

169 

wr  t 

6.5,"st.ds.  test  sec.,C(38J 

6 

",C|  381  ,  " 

",D|  381 

1  70 

wr  t 

6,"  Keyboard 

Inputs" 

171 

wr  t 

6,"  - 

172 

f  mt 

5, 1 12 . 6 

173 

wr  t 

6 . 5 , "Atmospher ic  Pres  sure ,C | 39 ) 

1  n  .  Hq 

"  ,D1  39  1 

174 

wr  t 

6. 5, "Water  flow  rate  thru  #1, 

C(40| 

Ibm/hr 

"  ,  D  :  4  0  1 

175 

wr  t 

6. 5, "Water  flow  rate  thru  #2, 

C|41| 

Ibm/hr 

" ,  L ;  4 1 1 

176 

wr  t 

6. 5, "Water  flow  rate  thru  #3, 

C(421 

Ibm/ hr 

”,01421 

1  7 

wr  t 

6. 5, "Water  flow  rate  thru  #4, 

C|43| 

Ibm/hr 

"  ,  P  1  4  i  1 

1  78 

wr  t 

5  ^ «  - - - 

179 

wtb 

6,12 

Fig.  B  -  2. 


Continued 


Operator : 


1  bC  : 
IBl  : 

j-  6  i  : 

18  i  : 
184; 
18S: 
186  : 
187  : 
168  : 

189  : 

190  : 

191  : 
192; 
19J  : 
194  ; 
195: 

196  ; 

197  : 

198  : 

1  99: 
200  : 
20  1  : 
20  2  ; 
20  i  : 
204  ; 
20  5: 
206  : 
20  7  : 
20  8: 
20  9  : 
210: 
^11: 
212: 
213: 
214; 

2  15: 
2  16: 
2  17: 
218: 
B  1  9  ; 
220  : 
221  : 
222: 
22  3: 
224  ; 
2  2  5  : 
2  26  : 
2  2  7  - 
2  2  8  : 
229  : 
2  iO  ; 
^11: 
^32; 

2  M  : 
2  14  : 
2  15: 
2  16  ; 
2  ’  7  : 
2  38  : 
2  19: 


wf  t 
wr  t 
w  r  t 
wr  t 
wr  t 
w  r  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
w  r  t 
wr  t 
wr  t 
wr  t 
f  mt 
wr  t 
wt  t 
wr  t 
wr  t 
sfg 
w  r  t 
wr  t 
wt  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wt  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
wr  t 
w  t  t 


6, 
6 , 
6 , 

6  , 


Date; 

Type  of  Put :  " , XS 


,AS 


Hednced  Run  Data 


6.5, 

"Air  density  at  orifice 

i  r 

t  ' 

6.5, 

"Air  viscosity  at  orifice 

i  i  iT 

■  '  r.  t  - 

I 

6.5, 

"Air  mass  flow  rate 

U'.T 

n  t 

6.5, 

"Average  inlet  air  temperature 

deu 

.  b 

6.5, 

"Average  outlet  air  temperature 

iJ  eg 

6.5, 

"Air  temper  a ure  drop 

d  r-Q 

.  r' 

6.5, 

"Specific  heat  of  air 

BTU, itT 

-deg 

.  [  ■' 

6.5, 

"Heat  transfer  from  air 

BTU  / 

h  I 

6.5, 

"Heat  transfer  to  water 

thru  #4 

B  Tv  - 

h  t 

6.5, 

"Heat  transfer  to  water 

thru  »  3 

HTU  •' 

hr 

6.5, 

"Heat  transfer  to  water 

thru  #2 

P.  TU  / 

h  t 

6.5, 

"Heat  transfer  to  water 

thru  1 1 

BTi:  / 

hr 

6.5, 

"Total  Heat  transfer  to 

water 

BTi) /.hr 

5  ,  f  1  2 . 6 

6.5, 

"Upstream  temp. drop  across 

tef Ion 

deg 

.  i' 

6.5, 

"Axial  conduction  from 

Ups 

t  rcan' 

BTL'/ 

h  r 

6.5, 

"downstream  temp. drop  across  teflon 

deg 

.  F 

6.5, 

"Axial  conduction  from  < 

downs  tream 

BTU/ 

hr 

14 

6.5, 

"Energy  Balance 

6.5, 

"Outlet  Air  Temperature 

at 

#1 

6f?q 

.  F 

6.5, 

"Heat  Tran.  Coeff.-l 

BTU/sq . f t-hr 

-deg 

.  F 

6.5, 

"Nusselt  Number-1 

6.5, 

"Prandtl  Number-1 

6.5, 

"Reynolds  Number-1 

6.5, 

"Outlet  Air  Temperature 

at 

#2 

deg 

.  F 

6.5, 

"Heat  Tran.  Coeff.-2 

BTU/sq. ft-hr 

-deg 

.  F 

6.5, 

"Nusselt  Number-2 

6.5, 

"Prandtl  Number-2 

6.5, 

"Reynolds  Number-2 

6.5, 

"Outlet  Alt  Temperature 

at 

#3 

dea 

.  F 

6.5, 

"Heat  Tran.  Coeff.-3 

BTU/sq. ft-hr 

-deg 

.F 

6.5, 

"Nusselt  Number-3 

6.5, 

"Prandtl  Numbec-3 

6.5, 

"Reynolds  Number-3 

6.5, 

"Outlet  Air  Temperature 

at 

#4 

deg 

.  F 

6.5, 

"Heat  Tran.  Coef f . -4 

BTU/sq .ft-hr 

-deg 

.  F 

6.5, 

"Nu.sselt  Number-4 

6.5, 

"Prandtl  Number-4 

6.5, 
6  .  "  - 

"Reynolds  Number-4 

boeptdsp  "Insert  tape  and  cont i nue " ; stp 
ent  "Track  1  or  0?",T;trk  T 
ent  "New  tape?  Yes= 1 : Ng=0 "  ,N 
if  N  =  1  ;  q  t  o  2  34 


ent  !’ 

fdf  P;idf  0,C',0,C 
if  O=0:qto  234 
Ptl»F;qto  231 
firt  6,f2.Q,fi.O 

wrt  6. 6, "Tape  Track=",T,"  Tape  File=",P 

asp  "Ready  to  mark  and  tape";stp 

if  N  =  1  ;  ] m p  2 

fdf  P 

mr  k  1,936 


"  ,  X  4 
,  X  . 

"  ,  X  6  ' 

"  ,  X  ;  , 

" ,X, a  ' 

"  ,  X  ,  1  0  , 

" ,  X :  i  2 
"  ,  X  ,  I  4  i 
" , X :  16; 

"  ,  X  ;  1  7  j 

" , Dl 4" 1 -D 1 46 i 
"  ,  X  ,  9  i 

" ,Bl 44 I-D| 451 
" ,  X ;  1 1  i 

"  ,  X  :  1  7  J  /  X 1  8  i 
",  A  i  1  I 
",  A  [  6  1 
"  ,A  iS  I 
"  ,A  i  10  1 
"  ,A|  11  ) 

",11  in 

",  B  I  6  I 
",Bf8| 

"  ,B(10  1 
"  ,B|  11  1 
"  ,  F  f  1  I 
"  ,  F  [  6  1 
"  ,F(8I 
"  ,F(10  1 
"  ,F|  11  ) 

"  ,  G  ( 1  1 
",  G  I  6  1 
"  ,G|8  1 
"  ,G110  1 
" ,0111] 


B  -  2.  Continued 


J  -i  (  :  ref  P  ,  A  S  ,  B  S  ,  X  S  ,  X  :  *  i  ,  D  ;  *  i  ,  A  (  *  j  ,  b  (  *  ;  ,  r  .  ♦  ; 
J41:  dsp  "Tape  recorded wth  6, 12;^ to  it 
14.:  end 
2  4  i  ;  "  n  1  r  m  "  : 

2  4  4  :  p i - p  3 

t4S;  .  e.'il  +  i,  42p-'i|  (plepj!)/2  +  459.67)-2.93e-9 
24d:  if  apE  I  p  3- (  pi -P‘3,/ p4 /p2)  )  <  .  1  ;q  t  o  24b 

24  7:  pl-p5/  p4  ,  p2-p 3;q to  2  45 

2  48 :  p8* p5*  In (  (pi-p7) /(pl-p6)  )/. 990256 (p3-p7 

249;  1/  ( 1,  p8-.004  78  3  )*p9 

250:  1. 3466 e -2+2. 2323e-5(pl+p3)/2»pi0 

251:  2.t7p9/12,'plO-pll 

2  5/ :  2 . 628  6e- 3 ( ( pl+p3) /2  +  4  59 . 67 ) " 1 . 5-  (  (p 1  tp3 

25  3;  pl2*p4/pl0->pl  3 

254:  5  .  72  24  2  5p2,  pi  2 -I'l  4 

255:  ret 


Fig 


B  -  2.  Concluded. 
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Table  C  -  1.  Heat  transfer  data  for  smooth  tube. 
TYPE  OF  INSERT:  No  Insert 


Table  C 


Heat  transfer  data  for  Type  A  insert  in  vertical  and 

horizontal  orientations. 


TYPE  OF 

INSERT: 

Type  A  in 

vertical  orientation 

Pr 

NUp 

Pr 

Rep 

NUd 

5916 

57.90 

0.706 

13567 

105.62 

0.707 

6228 

70.48 

0.707 

14050 

113.73 

0.707 

6515 

69.82 

0.707 

14955 

110.46 

0.705 

6752 

67.00 

0.707 

15561 

120.82 

0.706 

7049 

64.76 

0.706 

16137 

117.82 

0.707 

7408 

76.21 

0.707 

16684 

123.53 

0.707 

7735 

75,69 

0.707 

17619 

125.53 

0.705 

8008 

71.63 

0.707 

17809 

121.62 

0.705 

8105 

69.37 

0.705 

18318 

136.23 

0.706 

8488 

82.55 

0.707 

18524 

139.24 

0.706 

8849 

81.87 

0.707 

18977 

132.84 

0.707 

9039 

76.22 

0.705 

19203 

134.09 

0.707 

9159 

77.86 

0.707 

19615 

147.12 

0.707 

9450 

89.62 

0.707 

19809 

127.72 

0.707 

9840 

88.49 

0.707 

20079 

135.91 

0.705 

10180 

83.96 

0.707 

20809 

146.12 

0.706 

10244 

81.82 

0.705 

21509 

143.05 

0.707 

10287 

85.59 

0.705 

22192 

153.03 

0.707 

10720 

97.28 

0.706 

22446 

148.43 

0.705 

10724 

95.88 

0.707 

23267 

159.75 

0.706 

11143 

96.27 

0.707 

23959 

153.26 

0.705 

11170 

95.45 

0.707 

24044 

153.56 

0.707 

11543 

103.58 

0.707 

24786 

162.92 

0.706 

11558 

90.32 

0.707 

24793 

165.18 

0.707 

12533 

97.36 

0.705 

25570 

156.82 

0.707 

13064 

107.08 

0,706 

26332 

168.12 

0.707 

TYPE  OF 

INSERT: 

Type  A  in 

horizontal 

orientation 

Rep 

Pr 

Reo 

Nuq 

Pr 

5876 

_ 

55.53 

0.706 

13866 

106.74 

0.707 

6190 

67.86 

0.707 

14349 

113.89 

0.707 

6480 

66.24 

0.707 

15258 

109.87 

0.705 

67  34 

73.01 

0.707 

15876 

125.66 

0.706 

6813 

62.03 

0.705 

16464 

118.61 

0.707 

7187 

75.96 

0.707 

17018 

130.09 

0.707 

7532 

74.95 

0.707 

17157 

111.36 

0.705 

783  3 

82.11 

0.707 

17808 

132.81 

0.706 

8836 

74.06 

0.705 

18446 

127.04 

0.707 

9291 

89.72 

0.707 

18871 

127.53 

0.705 

9715 

88.0  3 

0.707 

19051 

136.95 

0.707 

10091 

95.20 

0.707 

19592 

147.66 

0.706 

12808 

9  3.45 

0.705 

20290 

141.51 

0.707 

1  3  346 

1  10.90 

0.706 

20954 

154.02 

0.707 
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Table  C  -  3.  Heat  transfer  data  for  Type  B  insert  in  vertical  and 
horizontal  orientations. 


TYPE  OF 

INSERT: 

Type  B  in 

vertical  orientation 

Reo 

NUd 

Pr 

ROd 

NUd 

Pr 

6058 

55.44 

0.705 

13368 

91.47 

0.707 

6359 

58.14 

0.707 

13805 

96.41 

0.707 

6626 

61.08 

0.707 

14901 

10201 

0.705 

6870 

65.20 

0.707 

15447 

98.30 

0.706 

7059 

62.35 

0.705 

15957 

103.27 

0.707 

7405 

65.54 

0.706 

16473 

110.45 

0.707 

7719 

68.75 

0.707 

17510 

115.80 

0.705 

8008 

72.81 

0.707 

18128 

108.84 

0.706 

8016 

63.94 

0.705 

18702 

115.17 

0.707 

8378 

64.79 

0.706 

19290 

125.76 

0.707 

8709 

70.22 

0.707 

19982 

127.34 

0.705 

9022 

73.73 

0.707 

20542 

130.38 

0.706 

9126 

69.30 

0.706 

20648 

118.26 

0.706 

9480 

71.30 

0.707 

21192 

125.67 

0.706 

9810 

74.48 

0.707 

21266 

124.81 

0.707 

10133 

81.62 

0.707 

21793 

127.78 

0.707 

10137 

7Z66 

0.705 

21899 

135.48 

0.707 

10312 

77.86 

0.705 

22391 

138.33 

0.707 

10554 

74.40 

0.707 

22878 

141.52 

0.705 

10698 

76.62 

0,706 

23596 

131.84 

0.706 

10933 

77.72 

0.707 

24263 

147.61 

0.705 

11064 

82.04 

0.707 

24266 

138.89 

0.707 

11290 

8Z90 

0.707 

24952 

149.57 

0.707 

11438 

88.99 

0.707 

24992 

136.37 

0.706 

12474 

88.01 

0.705 

25671 

142.58 

0.707 

12934 

85.10 

0.706 

26366 

152.38 

0.707 

TYPE  OF 

INSERT: 

Type  B  in 

horizontal 

orientation 

RCd 

NUd 

Pr 

NUd 

Pr 

5635 

43.83 

0.705 

12589 

84.37 

0.705 

5898 

51.4  3 

0.706 

13108 

94.09 

0.706 

6147 

50.% 

0.707 

13593 

90.25 

0.707 

6374 

54.65 

0.707 

14040 

94.78 

0.707 

7025 

^7.73 

0.705 

14317 

92.94 

0.705 

7325 

56.77 

0.706 

14888 

10232 

0.706 

7616 

56.50 

0.707 

15422 

98.17 

0.707 

7887 

60.05 

0.707 

15922 

103.13 

0.707 

8582 

66.37 

0.705 

16832 

105.89 

0.705 

8986 

74.28 

0.706 

17481 

115.74 

0.706 

9359 

73.26 

0.707 

18091 

111.67 

0.707 

9697 

76.53 

0.707 

18667 

117.91 

0.707 

10391 

71.06 

0.705 

18778 

113.58 

0.705 

10840 

80.04 

0.706 

19472 

123.89 

0.706 

11256 

77.85 

0.707 

20129 

120.13 

0.707 

11637 

82.72 

0.707 

20752 

126.16 

0.707 

o! 

OJ 

05  1 
1 
1 

_ 

Pr 

1 

1 

1  -v 

(D 

lO 

1 

1 

1 

Pr 

5797 

40.17 

0.706 

13028 

_ _ 

64.58 

0.706 

6003 

36.90 

0.707 

13390 

65.37 

0.707 

6188 

39.28 

0.707 

14846 

77.87 

0.704 

6370 

40.78 

0.707 

15304 

69.57 

0.706 

6905 

48.27 

0.705 

15718 

72.65 

0.706 

7184 

44.22 

0.706 

16132 

74.54 

0.707 

7428 

45.98 

0.707 

17457 

87.36 

0.704 

7660 

46.73 

0.707 

17957 

78.80 

0.706 

7747 

50.78 

0.705 

18417 

81.40 

0.706 

8028 

46.73 

0.706 

18882 

84.16 

0.707 

8281 

48.89 

0.707 

20026 

94.44 

0.705 

8532 

51.08 

0.707 

20549 

85.64 

0.706 

8893 

53.29 

0.705 

21039 

89.32 

0.706 

9200 

48.88 

0.706 

21541 

91.06 

0.707 

9476 

50.43 

0.707 

22427 

10Z18 

0.705 

9714 

58.35 

0.704 

22626 

102.41 

0.705 

9747 

51.98 

0.707 

22987 

92.56 

0.706 

10014 

58.20 

0.706 

23182 

93.24 

0.706 

10048 

5Z22 

0.706 

23507 

94.05 

0.706 

10322 

54.24 

0.706 

23701 

95.58 

0.706 

10349 

55.14 

0.706 

24034 

96.46 

0.707 

10603 

56.08 

0.707 

24228 

97.01 

0.707 

10650 

56.27 

0.707 

25441 

110.87 

0.705 

10879 

56.98 

0.707 

26023 

101.19 

0.706 

12261 

68.93 

0.704 

26573 

103.69 

0.706 

12664 

61.55 

0.706 

27140 

107.75 

0.707 

table  C  -  5.  Friction  data  for  heated  runs. 


TYPE  CF  INSERT:  No  Insert  TYPE  OF  INSERT:  Type 


Rep 

f 

f 

5381 

0.0104 

6353 

0.0920 

5408 

0.0116 

7550 

0.0904 

6382 

0.0105 

8650 

0.0880 

6385 

0.0088 

%27 

0.0876 

7352 

0.0096 

10923 

0.0873 

7353 

0.0098 

10924 

0.0953 

8227 

0.0085 

13303 

0.0924 

10754 

0.0079 

15834 

0.0915 

14218 

0.0073 

18632 

0.0904 

17722 

0.0068 

18836 

0.0839 

19786 

0.0067 

21147 

0.0888 

23254 

0.0064 

23638 

0.0879 

26073 

0.0062 

25162 

0.0870 

Appendix  D:  Sample  Calculations  and 


Uncertainty  Analysis 


Given  below  are  equations  used  for  calculation  of  experimentally- 
determined  quantities  necessary  for  determination  of  the  surface  con¬ 
vection  coefficient  h  in  one  section  of  the  calorimeter.  Included  is 

a 

an  uncertainty  analysis  of  each  calculated  quantity  obtained  using  the 
method  of  Kline  and  McClintock  [13];  the  measurement  uncertainties  used 
in  this  method  are  based  primarily  on  instrumentation  manufacturers 
specifications  with  odds  of  20  to  1.  Measurement  uncertainties  are 
given  in  Table  Dl.  Use  of  necessary  conversion  factors  to  obtain  com¬ 
patible  units  is  assumed. 


Table  Dl.  Uncertainties  of  measured  quantities. 


Air  temperature 
Water  temperature 
Wall  temperature 


Water  flow  rate 


Properties  of  air 

(C  ,  k  ,  p  ,  P  ) 
p ,  a  a  a  r 

Pressure 


±0.4°F 


±0.4°F 


±0.75  Ibm/hr 


0.01  in.  H^O 


Calculation  of  air  density  for  orifice  meter 


Calculation  of  uncertainty  in  air  mass  flow  rate 


W-  = 
m 

a 


'Ap  p  W- 
or  a  C, 


v2  /C^  /Ap 


or 


2  /s 

W  1  + 

3  3  > 


The  result  is  W*  =  2.58  Ibm/hr 


Calculation  of  energy  gain  by  water 


q  =  mc  (T  -T  ) 
w  w  p,w  w2  wl 


The  inputs  are:  =  65.76  Ibm/hr 


c  =  1  BTU/lbtn-°F 
p,w 


(T  ,  -  T  ,)  =  22.67°F 
w2  wl 


The  result  is  =  1490.75  BTU/hr 


Calculation  of  uncertainty  in  energy  gain  by  water 


q^  w  p,w  AT^  p,w  w  m  w 


where  AT  refers  to  the  temperature  drop  (T  „  -  T  ) 
w  ^  w2  wr 


The  result  is  W  =57.25  BTU/hr 


Calculation  of  energy  gain  by  conduction  in  calorimeter 

/AT  A 


^c  =  ^ 


The  inputs  are 


t 

A 


1.77  BTU/hr-ft-°F 
2 


AT.  = 


t 

Ax 


0.05  ft 
92.78°F 
1.4  in . 


The  result  is  q  =  70.44  BTU/hr 
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Calculation  of  uncertainty  in  energy  gain  by  conduction 


Considering,  W^, 

The  result  is  W  =  0.75  BTU/hr 


Calculation  of  corrected  energy  gain  by  water 
\?,c  ^c 

The  inputs  are;  =  1490.75  BTU/hr 

=  70.44  BTU/hr 

The  result  is  q^  =  1420.31  BTU/hr 

Calculation  of  uncertainty  in  corrected  energy  gain 

W  =  W  +  W 
q.c  q^,  c 

=  58.0  BTU/hr 

Calculation  of  outlet  air  temperature 


c 

p,a 


m 

a 


The  inputs  are;  ^al  ~  318.52°F 

q  =  1420.31  BTU/hr 

w ,  c 

C  =  0.2473  BTU/lbm-°F 

p,a 

m  =  125.22  Ibm/hr 

a 


The  result  is  T 


=  272.66°F 


The  inputs  are: 


=  1420.31 


A  =1.02  ft 

c 

ATl^td  =  165.36°F 


The  result  is  U  =  8.76  BTU/hr-ft  -°F 

Calculation  of  uncertainty  in  overall  coefficient 


2  AT 


Calculation  of  average  heat  transfer  coefficient  of  air 

1 


h  = 
a 


1 

U  2k 


LMTDi 


The  inputs  are,  D  =  2.67  in. 


Dj.  =  2.80  in. 


=  25  BTU/hr-ft-°F 
U  =  8.76 


The  result  is  h^  =  8.78  BTU/hr-ft  -°F 


Calculation  of  uncertainty  in  average  heat  transfer  coefficient 


65 


The  result  is  W-  =  0.4015  BTU/hr- f t^-°F 
ha 


Calculation  of  Nusselt  Number 


Nup 


h  D 
a 

k 


a 


The  inputs  are; 


The  result  is  Nu^  = 


=  8.78  BTU/hr-ft^-°F 

D  =  1 . 67  in. 

k  =  0.0201 
a 

97.36 


Calculation  of  Reynolds  Number 
,  m 

Re-  -  —  — 

D  n  u  D 
a 

The  inputs  are:  m  =  125.22  Ibm/hr 

a 

p  =  0.05717  Ibm/ft-hr 
a 

D  =  2.67  in. 

The  result  is  Re^^  =  12,533 


Calc ulat ion  of  uncer tainty  in  Reynolds  number 


The  result  is  287 

Reo 


Calculation  of  correction  for  pressure  drop  due  to  temperature 
di f  ference 


A  -  a  (  1  1 

^^c  A  Ip  p 

X  \  al  a2 


where  subscripts  1  and  2  refer  to  inlet  and  exit  locations 

The  inputs  are:  m  =  125.22  Ibm/hr 

a 

p  ,  =  0.04536  Ibm/ft^ 
al 

p  .  =  0.06145  Ibm/ft^ 
a  2 

A  =  0.0389  ft^ 

X 

Ap  =0.03  in.  of  H„0 
c  z 


Calculation  of  the  uncertainty  in  correction 


57 


1 r ulation  of  uncertainty  in  corrected  pressure  drop 


w  -  w.  +  w. 

Ap^  Ap  Ape 


The  result  is  W  =  0.015  in.  HO 

Ap^  2 


Calculation  of  friction  factor 


a  V 


D 

2  L 


2 

71 

16 


P  D" 
a 


The  inputs  are; 


Ap^  =  0.44  in.  H^O 

m  =  125.22  Ibm/hr 
a 

p  =  0.05654  Ibm/ft' 
a 


D  =  2.67  in. 

L  =  7 1 .  in . 

The  result  is  f  =  0.09243 

Cal culation  of  uncertainty  in  friction  factor 


APf^a 


-Ap  3  D  W  '  2 

+  1  ^  - ±  I  + 

•2  2 
m  L 
a 


,  2  1/2 

-2  p-p  D  w* 
f^a  m 


The  result  is  = 


0.0052 


LIST  OF  SYMBOLS 


A 

A 


p,a 


p,w 


t 

fn. 


fn. 


area  for  conduction  through  teflon  spacer 
empty  tube  surface  area  for  calorimeter  section 
cross  sectional  flow  area  of  calorimeter  section 
coefficient  of  discharge 

specific  heat  at  constant  pressure  of  air 
specific  heat  at  constant  pressure  of  water 
inside  diameter  of  calorimeter  pipe 

diameter  of  location  of  thermocouples  in  calorimeter  pipe  wall 
function  of  D  and  H  in  Eq.  (1)  from  [8] 
function  of  D  and  H  in  Eq.  (2)  from  [8] 


H 

axial  length  for  one  full  twist  in  Eqs . 

(1)  and 

1 

h 

a 

average  heat  transfer  coefficient  for  ai 

r 

k 

a 

thermal  conductivity  of  air 

k 

s 

thermal  conductivity  of  calorimeter  wall 

1 

thermal  conductivity  of  teflon 

L 

overall  length  of  calorimeter  (distance 

between 

L 

c 

heat  transfer  length  of  a  calorimeter  se 

ction 

♦ 

m 

a 

mass  flow  rate  of  air 

ni 

w 

mass  flow  rate  of  water 

P 

pressure 

< 

^amb 

ambient  air  pressure 

Psl 

static  air  pressure 

Ap  measured  pressure  drop  over  calorimeter  length 


Ap^  correction  to  measured  pressure  drop  over  calorimeter  length 
Ap^  corrected  pressure  drop 
Ap^^  pressure  drop  across  orifice 
energy  given  up  by  air 

q^  energy  loss  by  conduction  in  talorimcLer 

q  energy  gain  by  water 

w 

q  corrected  energy  gain  by  water 

w ,  c 

K  ideal  gas  constant 

inlet  air  temperature  to  a  calorimeter  section 

Ta2  ‘iif'  temperature  from  a  calorimeter  section 

average  air  hulk  temperature 

average  wall  (surface)  temperature  in  a  calorimeter  section 
static  temperature  at  orifice  meter 
j  inlet  water  temperature  to  a  calorimeter  section 
^s2  water  temperature  from  a  calorimeter  section 

AT^  temperature  difference  across  teflon  spacer 
U  overall  heat  transfer  coefficient 

AX  teflon  spacer  thickness 

W  uncertainty,  subscripted  for  each  quantity  for  which  the 
uncertainty  is  known  or  is  to  be  calculated 

Dimensionless  parameters 

f  friction  factor,  defined  in  F,q.  (5) 

Gr  Grashof  number 

NUp  Nusselt  number  based  on  inside  pipe  diameter 

Pr  Prandtl  number 

Reynolds  number  based  on  inside  pipe  diameter 
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